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g : ABSTRACT 
O 

We present rapid-cadence Transition Region And Coronal Explorer (TRACE) 
observations which show evidence of a filament eruption from active region 
00 ' NOAA 10696, accompanied by an X2.5 flare, on 2004 November 10. The eruptive 

filament, which manifests as a fast coronal mass ejection some minutes later, rises 
as a kinking structure with an apparently exponential growth of height within 
TRACERS field of view. We compare the characteristics of this filament eruption 
^ ■ with MHD numerical simulations of a kink-unstable magnetic flux rope, find- 

■ ing excellent qualitative agreement. We suggest that, while tether- weakening by 
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breakout-like quadrupolar reconnection may be the release mechanism for the 



O ■ previously confined flux rope, the driver of the expansion is most likely the MHD 

• helical kink instability. 
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1. Introduction 

Elucidating the release and driver mechanisms of both coronal mass ejections (CMEs) 
and filament eruptions is a prominent issue in current solar physics research. Several physical 
mechanisms are in competition. In the following, we briefly describe three of them, selected 
because the data we analyze are able to impose constraints upon them. 

The anchorage of magnetic field in the photosphere (and below) provides a strong stabi- 
lizing mechanism against the magnetic stresses which build magnetic pressure. For example, 
a stressed core field receives a strong stabihzing force from the magnetic tension of its over- 
lying magnetic arcade, whose field lines are analogous to tethers. Cutting these tethers, via 
magnetic reconnection, is one possible way of allowing the core field to erupt. Starting with 
a simple, bipolar, sheared magnetic arcade, van Ballegooijen & Martens (1989) proposed 
that converging footpoint motions will force magnetic reconnection between the arcade field 
lines, progressively transforming the arcade core into a twisted flux tube. Such reconnection 
weakens the photospheric anchorage and, at some point, the remaining arcade tethers arc 
expected to be too weak to prevent the eruption of the stressed core. In this approach, mag- 
netic reconnection is envisioned as working at low altitude, below the stressed fleld which 
will later erupt. This approach has been developed, and compared with observations, by 
Moore & Roumeliotis (1992). 

An alternative to this approach is one where reconnection above the stressed fleld is 
responsible for the removal of overlying tethers, namely the magnetic "breakout" model 
(Antiochos, DeVore, & Klimchuk 1999; Lynch et al. 2004). This requires a large-scale, 
quadrupolar magnetic field configuration. The core of the central arcade is assumed to be 
increasingly sheared by photospheric motions parallel to the magnetic polarity inversion line 
(IL). This core then expands with a nearly horizontal current sheet being formed around 
the X-point of the quadrupolar configuration. Magnetic reconnection, at this current sheet, 
progressively removes the field lines above the central, sheared arcade by transferring these 
magnetic connections to the lateral arcades (the two side-lobes of the quadrupolar configu- 
ration). At some point in this evolution, the feedback between the core expansion and the 
removal of the field fines above becomes explosive and the eruption of the core field follows. 
Gary & Moore (2004) analyzed an eruptive fiare which has all the main characteristics of 
the breakout model. A lateral version of this breakout model is also possible: here the 
most stressed field is in one of the side-lobes of the quadrupolar arrangement. In this case, 
reconnection only weakens the stabilizing upper arcade (by transferring it to larger scales). 
Observed flares demonstrating the characteristics of this mechanism have been analyzed by 
Aulanier et al. (2000), Gary & Moore (2004) and Sterling & Moore (2004). 

Another possible eruption mechanism is based on the helical kink instability (KI) of a 
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coronal magnetic flux rope. The KI occurs if the twist, a measure of the number of windings 
of the field fines about tfie rope axis, exceeds a critical value (Hood & Priest 1981), leading 
to a fielical deformation of tfie flux rope axis. In coronal applications of curved and line-tied 
flux-rope models, tfiis deformation results in an exponentially growing fieigfit of tfie flux 
rope's apex in tfie linear pfiase of tfie instability, as sfiown in numerical simulations (Fan 
& Gibson 2003; Torok, Kliem, & Titov 2004). Based on tfiese simulations, Torok & Kliem 
(2005) and Fan (2005) succeeded in modelling tfie ejection of a coronal flux rope from the 
Sun. Rust & LaBonte (2005) also presented observational evidence that the KI is the driver 
in a number of filament eruptions. 

Are tfiese CME models mutually exclusive, tfien, or is a combination at work on tfie 
Sun? In tfiis letter, we investigate a well observed event tfiat occurred on 2004 November 10 
in order to look for clues. 

2. Observations &c Reduction 

We differentially rotated 162 full-disk magnctograms, taken witfi tfie Solar and Hclio- 
spheric Observatory's Micfielson Doppler Imager (SoHO/MDI; Scfierrer et al. 1995) between 
2004 November 01 and 2004 November 10, to tfie central meridian, and analyzed subflelds 
centered around Region NOAA 10696. Tfiis allows tfie deduction of a negative sign of mag- 
netic fieficity for tfie flux emerging in tfiis region (Lopez Fuentes et al. 2003). During tfie 
evolution of tfiis region, tfiere is continuous input of magnetic energy and fieficity, wfiicfi 
renders it prolific in flares and CMEs. Here, we fiave investigated tfie eruption of a filament 
in tfie centre of Region 10696 on 2004 November 10 (Fig. 1), wfiicfi seems to be part of tfiis 
longer evolution and led to a very fast CME. We note tfiat, in tfie vicinity of tfie filament 
(fiereafter indicated by /), tfiere is a westward sfiearing motion of large positive-polarity 
elements, parallel to tfie IL. 

Tfie Transition Region And Coronal Explorer {TRACE; Handy et al. 1999) observed 
tfiis eruption primarily in tfie 1600 A band. From 01:21:56 UT, tfiese observations were 
made at a cadence of ~34 s until 01:57:30 UT, and tfiereafter at a mean cadence of 3.7 s. In 
order to sfiow tfie pre-eruption filament, we also present an image (Fig. lb) from tfie 195 A 
filter of SoifO's Extreme-ultraviolet Imaging Telescope (EIT; Delaboudiniere et al. 1995) at 
01:12:50 UT. 

For consistency, all data sfiown in Figs, la to le (and tfie movie ktrace.mpg) were 
differentially rotated back to tfie time of tfie last MDI image before eruption at 01:39:02 UT. 

At approximately 01:51:40 UT, a fiigfily localized brigfitening occurred close to tfie IL 
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at (Fig. la; see corresponding lightcurve in Fig. 2). This brightening is approximately 
co-spatial with the emergence of a small bipole (hereafter indicated by EF), north of the IL 
(Fig. le); we assume that its lack of appearance until 03:15 UT is due both to MDl's long 
cadence and to the small scale of the bipole emerging in a strong, positive- flux environment. 

The initial brightening at was followed by near-simultaneous brightenings at this 
and three other sites (marked B~ , and D~ in Fig. la; superscripts indicate magnetic 
polarity) between ~01:53 UT and 01:58 UT. These four, ribbon-like brightenings precede the 
later flare ribbons, the latter being detected as sharp rises in 1600-A flux from -02:02 UT 
(Fig. 2). 

The eastern leg of / began to become visible in emission at 02:01:26 UT, with the full 
fllament visible by 02:02:50 UT. The fllament's curved axis had already developed a slight 
cusp shape near its apex by the latter time, which proceeded to develop further (Figs. lc,d). 
The eruption of this fllament is associated with the development of J-shaped flare ribbons 
on both sides of the IL (Fig. Ic). 

3. Interpretation 

3.1. Eruption topology 8z kinematics 

The initial brightening at indicates that the emerging bipole EF (Sect. 2) reconnects 
with the pre-existing coronal field to give new connectivities (Fig. 3a). The newly reconnected 
field lines above the filament / are longer than the original ones, and so provide less downward 
stabilizing tension. Moreover, as with all emerging flux, the bipole EF carries magnetic 
energy and hclicity and a fraction is transferred to the large, reconnected loops, rendering 
them more stressed and thus able to expand further. 

The four observed precursor ribbons are evidence that a quadrupolar-like reconnection 
occurred in this eruption. This may seem surprising in this magnetic configuration, formed 
mainly by a positive and a negative magnetic region, so similar to a bipolar magnetic arcade. 
However, both positive and negative regions are divided into separate magnetic concentra- 
tions, and the IL is curved. This type of configuration is known to have a complex topology, 
not necessarily in the classical sense (i.e., with separatrices) but in the broader context of 
Quasi-Separatrix Layers (QSLs; Demoulin et al. 1997). QSLs are very thin layers where 
current sheets, and subsequently 3D magnetic reconnection, can develop. The equivalent 
2D magnetic topology is quadrupolar, as shown in Fig. 3. The emergence of bipole EF 
and - on a longer time scale - the observed shearing motions in the vicinity of / (Sect. 2) 
both contribute to an expansion of the magnetic configuration supporting /. This can trig- 
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ger quadrupolar reconnection, which weakens the tethers of the overlying arcade. The four 
ribbons, and plausibly also the release of /, are then interpreted naturally by the lateral 
breakout mechanism described in Sect. 1 and Fig. 3b. 

The erupting filament, however, shows evidence of a kinked shape as soon as it becomes 
observable by TRACE. Moreover, the associated eruptive ribbons are J-shaped (Fig. Ic), 
just as expected with the QSLs computed for twisted magnetic configurations (Demoulin, 
Priest, & Lonie 1996). The separation of the ribbons along the IL, as well as their reversed- J 
shape, are characteristic of a negative helicity, again supporting the conclusion derived from 
the magnetic evolution (Sect. 2). The main origin of the eruption is therefore most likely 
the accumulated magnetic stress, which progressively drives the configuration toward being 
kink-unstable. 

In order to quantify the rise behaviour of the kinking filament with time, the projected 
position of its apex was determined by selecting co-ordinates at the apparent apex using 
a mouse cursor. We define the projected apex height as the distance, in the plane of the 
sky, between the apex and the point directly between the extremities of the filament as seen 
at 195 A (Fig lb). The accuracy of these positions was tested by a series of independent 
measurements. The errors in these measurements are estimated as the quadrature sum of 
trace's FWHM (0.5") and the 1-pixel error on the screen during measurement, in each 
direction: this yields a lo" error in height of 1". 

We then fit these height-time data (Fig. 4) with several growth laws (constant acceler- 
ation, linearly increasing acceleration, power-law, exponential) and find that it is difficult to 
distinguish between these models from the limited data in TRACE'S field of view. The most 
satisfactory fit is to an exponential function: 



where {ho + Oq) = 1.5 x 10"^ ± 6 x 10^ km is the height of the filament at 02:02:50 UT {t=0), 
Oq = 6.7 X 10^ ±2.7 X 10^ km, and r is the e-folding time (or 'growth time') of the exponential 
acceleration term. In equation (1), r is unaltered by projection effects, and is found to be 
best fit by T = 95 ± 2 s. The power-law and constant-acceleration models were both found 
to fit the data comparably well in terms of residuals, but these yielded start times near 
02:04:15 UT; both the TRACE and GOES observations suggest, however, that the eruptive 
process had already started by 02:02:50 UT. 



h{t) — hf) + aQ exp (t/r) 
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3.2. Numerical simulation 

The erupting filament / exhibits a clear helical deformation of its axis shape (Fig. 1) as 
well as an apparently exponential growth in apex height (Fig. 4). These observations suggest 
that the filament eruption itself is primarily driven by the KI of a magnetic fiux rope which 
contains the filament. 

To substantiate this assumption, we compare our observations with the numerical sim- 
ulation of Torok & Kliem (2005). They used a modified version of the coronal fiux rope 
equilibrium developed by Titov Sz Demoulin (1999) as the initial condition in 3D, ideal-MHD 
simulations. By prescribing a supercritical flux-rope twist and a sufficiently weak ambient 
magnetic fleld, the kink-unstable flux rope showed an initially exponential rise, followed by 
a rise with a high constant velocity. 

The Titov & Demoulin equilibrium can be regarded as a model of a bipolar active region 
with the flux rope supporting a fllament. The large-scale photospheric flux distribution is 
qualitatively similar to the distribution in Region 10696 before the eruption: two main 
polarity regions are divided by a curved IL (see Fig. 4 in Titov & Demouhn 1999). The flux 
rope orientation in the model is nearly aligned with the IL (Figs. If thru Ih), as is generally 
the case for filaments. 

In Fig. 1, we compare the evolution of the axis shape of the filament with the simulation. 
The point of view towards the model fiux rope is chosen such that: i) the IL is oriented as 
in the observation; and ii) the morphology of the filament and the fiux rope match as close 
as possible. The (left-handed) helical axis deformation is qualitatively very similar on both 
counts. The path of the flux-rope apex deviates from radial ascent by about n/A northward in 
projection, presumably in the direction of the weakest overlying field (as one can extrapolate 
from Fig. le, and as indicated in Fig. 3) 

By comparing the morphology at corresponding times, we also find that the observed 
helical deformation of the filament is initially stronger than in the simulation. This indicates 
that the magnetic field strength, at least in the coronal region just above the filament, 
decreases more slowly with height than in the model (see Figs. 1 and 5 in Torok & Kliem 
(2005) for a comparison of the helical fiux rope deformation in two configurations with 
different ambient fields) . Future quantitative studies may enable the estimation of the height 
profile of the overlying magnetic field strength. 

A scaling of the simulation data to the observed height-time profile, using an initial 
height of 15.5 Mm at 02:01:10 UT and an Alfven time of 11.8 s {va = 1300 km s"^), also 
yields excellent agreement (Fig. 4), further supporting the conclusion, from the fitting, that 
the rise was exponential in nature. 
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4. Conclusions 

The emerging bipole north of the IL, spatially coincident with a pre-emption brightening 
in TRACERS 1600 A filter, suggests an initiation role for tether-weakening in the release of 
an already highly twisted filament. However, this bipole is relatively small ($ ~ 4.0 x 
10^^ Mx) when compared with the flux of the active region ($ Ri 4.5 x 10^^ Mx), and its 
role in weakening the field above the filament's flux rope is likely to be rather limited. Its 
true contribution could rather be to bring a small amount of magnetic stress (or helicity), 
pushing the twisted configuration toward the critical value for eniption or triggering of the 
quadrupolar reconnection which leads to a strong tether weakening. Four bright ribbons 
are also observed before the filament eruption: they are the signatures of a quadrupolar 
reconnection taking place well above the filament, and so of a lateral breakout mechanism. 
In this for Hori et al. (2005), evidence of both tether- weakening and breakout is 

present in the eruption. Furthermore, the filament eruption itself appears to follow an 
exponentially growing form (cf. Gallagher et al. 2003; Shanmugaraju et al. 2003), suggestive 
of an instability. The apparent writhe of the filament, and hence of the host flux rope, points 
to the MHD helical kink instability as the likely driver behind the eruption. Moreover, the 
characteristics of the filament eruption have the same global characteristics as the MHD 
simulations of a kink-unstable magnetic flux rope. 

We conclude that the release of the accumulated magnetic stress in the corona can 
involve several mechanisms, three of them being present in the studied event. The magnetic 
tethers are weakened both at low heights and well above the filament. However, the main 
engine of the filament eruption itself has all the observable characteristics of a kink instability. 

The authors are very grateful to Mitch Berger and Jingxiu Wang for helpful discussions, 
and to Jason King for valuable help in data retrieval. Work supported by grants OTKA 
T038013 & OTKA T048961 (LVDG), DFG MA 1376/16-2 (BK) and computing time at 
NIC, Jiihch. 
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Fig. 1. — a) TRACE 1600 A image: overlaid are four boxes to show the areas of the 
quadripartite brightening before the impulsive phase of the flare, b) EIT high-resolution 
image at 195 A showing filament / in its pre-eruption position between (725", 115") and 
(755", 110"); white box shows the area covered by the adjacent inset. Inset shows filament 
traced by a solid white line, with the apparent footpoints of the filament marked by black 
crosses, c) and d) TRACE 1600 A images showing the evolution of the bright kinked filament. 
Insets as for (b) but with black hne tracing the kinked shape of the erupting structure. Note 
the reversed-J shape of the ribbons between the filament's apparent feet in c). e) The first 
MDI longitudinal magnetogram after the flare's impulsive phase. The black circle denotes 
a small patch of negative flux which has appeared north of the IL. f) thru h) magnetic 
field lines outlining the core of the kink- unstable flux rope at t=0, 21 and 25 Alfven times 
in the numerical simulation of Torok & Kliem (2005). The grayscale levels in the plane 
indicate positive (white) and negative (black) concentrations of the magnetic field strength 
component normal to the plane. Figs, a), c) and d) form part of ktrace .mpg. Fig. e) is part 
of kmdi .mpg, and Figs, f) thru h) are taken from ksimu.mpg in the electronic edition of the 
Astrophysical Journal. 



- 11 - 




01540 D1?I5 01;5O 01:5& OZflO 02i06 OilO S 
UT on 2EXH-NmmiHr-1Q 



Fig. 2. — Smooth dashed curve: GOES soft X-ray flux; scale given in GOES class on left 
vertical axis. Solid histogram curves: pixel- averaged TRACE 1600-A intensity lightcurves 
from the four boxed regions shown in Fig. la); scale on right vertical axis. The horizontal 
bar indicates the range of time covered by data in Fig. 4. 
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Fig. 3. — Main topological features of the magnetic configuration, with a simplified realiza- 
tion along the curved IL, folded into a 2D diagram. Dotted lines represent key separatrices; 
dashed lines, former connectivities later transformed by reconnection; continuous thick lines, 
new connectivities after reconnection; and continuous thin lines, continuing connectivites. a) 
Evolution in the initial phase of the eruption assuming a tether-weakening due to reconnec- 
tion with the new flux, EF . b) The eruption with reconnection both above (quadrupolar) and 
below the erupting filament. The bottom- right arrow in each panel indicates the strongest 
photospheric motions parallel to the IL. 
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Fig. 4. — Log-linear plot of the height-time profile of the filament apex in the plane of the 

sky. Vertical bars represent data points with la errors. The dashed line represents the best 
fit to the data by equation (1). The solid line represents the h{t) data of the simulation 
shown in Fig. 1. Lower panel shows residuals to the fit. 



